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Abstract 

In this paper we study the maximum number of hmit cycles that can 
bifurcate from a focus singular point pQ of an analytic, autonomous differen- 
tial system in the real plane under an analytic perturbation. We consider po 
being a focus singular point of the following three types: non-degenerate, de- 
generate without characteristic directions and nilpotent. In a neighborhood 
of po the differential system can always be brought, by means of a change 
to (generalized) polar coordinates {r,9), to an equation over a cylinder in 
which the singular point po corresponds to a limit cycle 79. This equation 
over the cylinder always has an inverse integrating factor which is smooth 
and non-flat in r in a neighborhood of 70. We define the notion of vanishing 
multiplicity of the inverse integrating factor over 79. This vanishing multi- 
plicity determines the maximum number of limit cycles that bifurcate from 
the singular point po in the non-degenerate case and a lower bound for the 
cyclicity otherwise. 

Moreover, we prove the existence of an inverse integrating factor in a 
neighborhood of many types of singular points, namely for the three types 
of focus considered in the previous paragraph and for any isolated singular 
point with at least one non-zero eigenvalue. 
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1 Introduction and statement of the results 

Let us consider a planar, real, analytic, autonomous differential system with a 
singular point which we assume to be at the origin, that is, we consider a differential 
system of the form: 

X = P{x,y), y = Q{x,y), (1) 
*The authors are partially supported by a DGICYT grant number MTM2005-06098-C02-02. 
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where P{x,y) and Q{x,y) are real analytic functions in a neighborhood U of the 
origin such that P(0,0) = (5(0,0) = 0. As usual, we associate to system ([T]) the 
vector field Xq = P{x,y)dx + Q{x,y)dy. We assume that the origin po = (0,0) 
is an isolated singular point, that is, there exists a neighborhood of it without 
any other singular point, and we assume that it is a monodromic singular point. 
Therefore, it is either a center (i.e. it has a neighborhood filled with periodic orbits) 
or a focus (i.e. it has a neighborhood where all the orbits spiral in forward or in 
backward time to the origin). 

We consider an analytic perturbation of system ([1]) of the form: 

X = P{x,y) + P{x,y,e), y = Q{x,y) + Q{x,y,e), (2) 

where e G is the perturbation parameter, < ||£:|| << 1 and the functions 
P{x,y,e) and Q{x,y,e) are analytic for {x,y) G U, analytic in a neighborhood of 
e = and P{x,y,0) = Q{x,y,0) = 0. We associate to this perturbed system (j2]) 
the vector field = {P{x,y) + P{x,y,e))dx + {Q{x,y) + Q{x,y,e))dy. 

We say that a limit cycle 7^ of system ([2]) bifurcates from the origin if it tends 
to the origin (in the Hausdorff distance) as e — 0. We are interested in giving 
a sharp upper bound for the number of limit cycles which can bifurcate from the 
origin po of system ([1]) under any analytic perturbation with a finite number p of 
parameters. The word sharp means that there exists a system of the form ([2]) with 
exactly that number of limit cycles bifurcating from the origin, that is, the upper 
bound is reahzable. This sharp upper bound is called the cyclicity of the origin po 
of system ([T]) and will be denoted by Cycl(A'e,po) all along this paper. 

As defined in [U |2T1 [29] , a Hopf bifurcation, also denoted by Poincare- Andronov- 
Hopf bifurcation, is a bifurcation in a neighborhood of a singular point like the 
origin of system ([1]). If the stability type of this point changes when subjected to 
perturbations, then this change is usually accompanied with either the appearance 
or disappearance of a small amplitude periodic orbit encircling the equilibrium 
point. 

We remark that in all the examples considered in this paper we take a perturbed 
system ([2]) with the origin as singular point, that is, -P(0,0,£:) = Q{0,0,e) = 0. 
Moreover, all the limit cycles 7e that bifurcate from the origin in our examples 
surround the origin. 

We consider systems of the form ([1]) where the origin po is a focus singular point 
of the following three types: non-degenerate, degenerate without characteristic di- 
rections and nilpotent (the definitions are stated below). The study of the cyclicity 
of a degenerate focus has been tackled in very few sources; we mention the papers 
[U 131 [ini [121 HSl 1221 [21] where some (partial) results can be found. In relation with 
normal forms and integrability of degenerate singular points, we cite jHl 1221 [25] . 
Our results do not establish that the cyclicity of this type of singular points is finite 
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but give an effective procedure to study it. In the three mentioned types of focus 
points, we will consider a change to (generalized) polar coordinates which embed 
the neighborhood U of the origin into a cylinder C = {(r, ^) G M x 5^ : |r| < 6} 
for a certain sufficiently small value of 5 > 0. This change to polar coordinates is a 
diffeomorphism in W — {(0, 0)} and transforms the origin of coordinates to the circle 
of equation r = 0. In fact, the neighborhood U is transformed into the half-cylinder 
in which r > 0, but we can consider the extension to the values in which r < 
by using several symmetries of the considered (generalized) polar coordinates. In 
these new coordinates, system ([1]) can be seen as a differential equation over the 
cylinder C of the form: 

To = ^c--*)- 

where J-'{r, 6) is an analytic function in C. The circle r = needs to be a particular 
periodic solution of the equation ([3]) and, therefore, ^-"(0, 9) = for all 9 E S^. 

Throughout the rest of the paper, we consider an inverse integrating factor 
V{r,9) of equation We recall that an inverse integrating factor of equation 
([3D is a function \^ : C — * M of class C^{C), which is non locally null and which 
satisfies the following partial differential equation: 

o9 or or 

We remark that since V{r, 9) is a function defined over the cylinder C it needs to 
be T-periodic in 9, where T is the minimal positive period of the variable 9, that 
is, we consider the circle 5^ = M/[0,T]. The function V{r,9) is smooth (C°°) and 
non-fiat in r in a neighborhood of r = 0. The existence of an inverse integrating 
factor V{r, 9) with this regularity is proved in using the result of [30], see also 
Lemma [2T] of the present paper. A characterization of when V{r,9) is analytic in 
a neighborhood of r = is given in [13] . 

Let us consider the Taylor expansion of the function V{r^9) around r = 0: 
V{r,9) = Vm{9)r'^ + 0{r"^~^^), where Vm{9) ^ for 9 E and m is an integer 
number with m > 0. As we will see in the following section, in fact, Vm{9) 7^ for 
all 9 E S^, cf. Lemma [201 We say that m is the vanishing multiplicity of V{r, 9) on 
r = 0. The aim of this paper is to show the correspondence between this vanishing 
multiplicity m and the cyclicity Cycl(A'e,po) of the origin pq of system ([T]). 

One of our hypothesis is that the origin of system ([1]) is a focus, and thus, we 
obtain that the circle r = is an isolated periodic orbit (i.e. a limit cycle) of 
equation ([3]). This hypothesis implies that there can only exist one inverse inte- 
grating factor V{r, 9) smooth and non-fiat in r in a neighborhood of r = 0, up to a 
nonzero multiplicative constant, see Lemma f2T\ The uniqueness of V{r,9) implies 
that the number m corresponding to the vanishing multiplicity of V{r, 9) on r = 
is well-defined. 
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The existence of V{r,6) gives the existence of an inverse integrating factor 
Vo{x,y) for system ([T]) by undoing the change to (generahzed) polar coordinates. 
We recall that Vq : W ^ M is said to be an inverse integrating factor of system 
([1]) if it is of class C^{U), it is not locally null and it satisfies the following partial 
differential equation: 

dVo{x,y) dVo{x,y) f dP{x,y) dQ{x,y) 

P{x, y) — 7^ + Q{x, y) = — + — 

ox oy \ ox oy 

The function Vo(x,?/), obtained from V{r,6) by undoing the change to polar co- 
ordinates, does not need to be smooth at the origin po since the change to polar 
coordinates is a diffeomorphism except in the origin. Besides, we are also inter- 
ested in the problem of the regularity of an inverse integrating factor Vo(x,?/) in a 
neighborhood of the origin of system ([T]), whenever it exists, and we also analyze 
this question, see Theorem [T6l and Corollaries [T2l and [T9l 

The zero-set of Vo(x, y), which we denote by V^~^(0) := {p ElA : Vo(p) = 0}, 
is formed by orbits of system ([1]) and usually contains those orbits which deter- 
mine the dynamics of the system: singular points, limit cycles and graphics, see 
[3 [151 [13 [IS]- If there exists an inverse integrating factor in a neighborhood of a 
limit cycle, then it is contained in V^~^(0), as it has been proved in |T9]. Since the 
origin po = (0, 0) of ([1]) is a focus, we have that Vo(0, 0) = 0, as proved in [71 [T7] . 
Moreover, the set Iq~^(0) contains all the singular points of ([1]) having at least one 
parabolic or elliptic sector in its domain of definition, as it is also proved in [3 [T7] . 

In [15], we proved that the vanishing multiplicity of an analytic y(r, 9) on the 
limit cycle r = coincides with the multiplicity (and, thus, the cyclicity) of r = 
as an orbit of equation ([3]). The statements and the proofs given in [15] can be 
repeated verbatim with the weaker assumption that V{r, 9) is smooth and non-flat 
in r in a neighborhood of r = 0, see Corollary [21] in the following section. The 
multiplicity of the limit cycle r = is related to the cyclicity of the origin of system 
([1]) by the change from polar coordinates. The aim of this work is to study the 
cyclicity of the origin of system ([1]) through the vanishing multiplicity of V{r, 9) 
on r = 0. We remark that the multiplicity m of r = as a limit cycle of equation 
([3D can also be established by successively solving the variational equations. How- 
ever, this method implies the computation of iterated integrals of non-elementary 
periodic functions. If we explicitly know an inverse integrating factor Vo{x,y) for 
system ([1]), then we have an inverse integrating factor V{r,9) of (|3]) and we can 
immediately know the value of m through the vanishing multiplicity of V(r, 9) in 
r = 0. On the other hand, since the existence of an inverse integrating factor 
V{r, 9), which is smooth and non-fiat in r in a neighborhood of r = 0, for equation 
([3]) is ensured (see Lemma [21]) . we have an alternative method to the variational 
equations to determine the value of m. 
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We are going to study the cyclicity of the origin of system ([T]) in the following 
three cases: the origin is a non-degenerate focus, the origin is a degenerate focus 
without characteristic directions and the origin is a nilpotent focus. The two first 
cases are transformed to an equation of the form ([3]) using polar coordinates and 
the latter case using the so-called generalized polar coordinates, see the definition 
below. We state our results for each type of coordinates, separately. The following 
Subsection 11.11 contains the results related to the non-degenerate and degenerate 
without characteristic directions case. Subsection 11.21 is devoted to the nilpotent 
focus. In Subsection 11.31 we analyze the existence (and regularity conditions) of an 
inverse integrating factor in a neighborhood of an isolated singular point of system 
([T]). We give some general results on an equation over a cylinder in Section [21 
Finally, Section [3] contains the proofs of our main results. 



1.1 A focus without characteristic directions 

We say that a focus at the origin of system ([T]) is non-degenerate if the linear part 
of system has complex eigenvalues of the form a ± f3i with a, /? G M and /3 7^ 0. 
After a linear change of coordinates and a rescaling of time, if necessary, system 
can be written in the form: 

X = (x - y + P2ix,y), y = X + Cv + Q2{x,y), (4) 

where C ^ 1^ and ^2(2;, y) and Q2{x, y) are analytic functions in a neighborhood of 
the origin without constant nor linear terms. 

We say that the origin of system ([T]) is a degenerate singular point if the deter- 
minant associated to the linear part of ([1]) is zero. We consider a system ([1]) of the 
form: 

X = Pd{x,y) + Pd+i{x,y), y = qd{x,y) + Qd+i{x,y), (5) 

where (i > 1 is an odd number, Pci{x,y) and qd{x,y) are homogeneous polynomials 
of degree d and Pd+i{x,y),Qd+i{x,y) G ?/)||'^+"^), that is, they are analytic 

functions in a neighborhood of the origin with order at least d + 1. We assume 
that Pd{x, y) + q^{x, y) ^ 0. When d > 1, the origin of system ([5]) is a degenerate 
singular point. 

A characteristic direction for the origin of system ([S]) is a linear factor in ]R[x, y] 
of the homogeneous polynomial xqd{x,y) — ypd{x,y). It is obvious that, unless 
xqd{x,y) — ypd{x,y) = 0, the number of characteristic directions for the origin of 
system ([5]) is less than or equal to d + 1. If there are no characteristic directions, 
then the origin is a monodromic singular point of system ([5]). We observe that 
the reciprocal is not true. A singular point with characteristic directions can be 
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monodromic. The origin of system ([T]) is a monodromic singular point if it is either 
a center or a focus, see for instance |H |T51 [21] for further information about mon- 
odromic singular points and characteristic directions. We assume that the origin 
of system ([5]) is a focus without characteristic directions. 

In the non-degenerate case {d = 1), that the cyclicity of a focus point is finite 
is well-known as well as several methods to determine it. The most usual method 
is to compute the first non-vanishing Liapunov constant and its order gives the 
cyclicity of the focus. Indeed, the same method allows to study the limit cycles 
which bifurcate from the origin in this case. When an inverse integrating factor is 
known, we give a shortcut in the study of the cyclicity as it can be given in terms 
of the vanishing multiplicity of the inverse integrating factor at the origin. 

As far as the authors know, the cychcity of a degenerate focus (rf > 1) of 
system ([5]) is not proved to be bounded. In fact, very few techniques appear in 
the literature to tackle the cyclicity of this type of singular points. Usually, polar 
coordinates are taken and the corresponding equation over the cylinder is studied. 
Our approach is to take profit from the knowledge of an inverse integrating factor 
to avoid the study of the differential equation over the cylinder. 

We remark that if d = 1 and the origin of system ([5]) is a focus point (without 
characteristic directions), then it can be written in the form (jl]). 

We use polar coordinates, {x,y) = (r cos 6, r sin^), in order to transform a 
neighborhood of the origin into the cylinder with period T = 2tt, and system 
into an ordinary differential equation of the form ([3]). 

In relation with system ([2]), an analytic perturbation field {P{x, y, e), Q{x, y, e)) 
is said to have subdegree s if {P{x,y,e),Q{x,y,e)) = 0{\\{x,y)\\'^). In this case, 
we denote by AfJ^' = {P{x, y) + -P(x, e))dx + (<5(a;, y) + Q{x, y, s))dy the vector 
field associated to such a perturbation. 

Theorem 1 We assume that the origin po of system ([5]) is monodromic and with- 
out characteristic directions. Let V{r, 6) be an inverse integrating factor of the 
corresponding equation which has a Laurent expansion in a neighborhood of 
r = of the form V{r, 6) = Vm{0) r™ + C(r'"+^), with Vm{0) ^0 andm eZ. 

(i) Ifm<Oormis even, then the origin of system is a center. 

(ii) // the origin of system is a focus, then m > 1, m is an odd number and 
the cyclicity Cycl(A'e,]9o) of the origin of system ([5]) satisfies Cycl(A'£,]9o) > 
(m + d) /2 — 1. In this casern is the vanishing multiplicity ofV{r,6) onr = 0. 

(ii.l) //, moreover, the focus is non-degenerate {d = 1), then the aforemen- 
tioned lower bound is sharp, that is, Cyc\{X^,Pq) = (m — l)/2. 
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(ii.2) // only perturbations whose suhdegree is greater than or equal to d are 
considered, then the maximum number of limit cycles which bifurcate 
from the origin is {m — l)/2, that is, Cycl(A'i'^',po) = {m — l)/2. 

The proof of this theorem is given in Section [3J 

Remark 2 As we will see in the proof of this theorem, if there exists an inverse 
integrating factor Vo{x,y) of system ([5]) such that Vo(rcos^, rsin^)/r'^ has a Lau- 
rent expansion in a neighborhood ofr = 0, then the exponents of the leading terms 
of Voir cos 6, rsin6)/r'^ and V{r,6) coincide. Thus, the vanishing multiplicity m 
can be computed without passing the system to polar coordinates. 

We provide several examples of application of Theorem [TJ 



Example 3 The following system 

X = -y {{2ij + l)x'^ + y^) + x^Xix^ + X2{x^ + y^)) , 

(6) 

y = a; (x' + (1 - 2/i)?/2) + a;2y(Aix2 + As (x^^ + 

where fi, Ai and A2 are real parameters, appears in |31j, where it is shown that 
the origin is a focus for a non semi-algebraic set of values of (/i, Ai, A2) and it is a 
center otherwise. 

We have that this system is written in the form ([5]) with d = 3 and it has no 
characteristic directions as xq^^x^y) — ypz{x,y) = {x^ + y'^y. Easy computations 
show that the function 

Vo{x,y) = e-'^+y'' (x + y y 

is an inverse integrating factor of system (jH]) which satisfies the hypothesis of our 
Theorem [T] with V{r,6) = e'"^^^"^ ^ r^. We remark that Vo(x, is not an analytic 
function in a neighborhood of r = unless /i = 0. As a consequence of Lemma [211 
which we prove in Section [3l we deduce that when the origin of system ([6]) is a focus, 
there are no analytic inverse integrating factors Vo(x, ?/) defined in a neighborhood 
of the origin. If it existed, its transformation to polar coordinates would produce 
an analytic inverse integrating factor V{r, 9) defined in a neighborhood of r = 
and different from ^(r, 6'), up to a multiplicative nonzero constant. In [31] it is 
shown that there exist values of (yU, Ai, A2), either with = or with /i 7^ 0, for 
which the origin of system ([6]) is a focus. Applying Theorem [H we deduce that 
whenever the origin of system ([B]) is a focus, then its cyclicity is greater than or 
equal to 2. In the particular case of a perturbation of subdegree greater or equal 
than 3, the maximum number of limit cycles that bifurcate from the origin is 1. 
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Example 4 Let k, s be integers such that s > 2A; > and consider the following 
differential system: 

X = -y{x^ + yY + xRs{x,y), y = x {x^ + y^ f + y Rs{x,y), (7) 

where Rs{x, y) is a homogeneous polynomial of degree s. The origin of this system 
is a monodromic singular point since there are no characteristic directions. We 
take polar coordinates x = r cos 9, y = r sin ^ and system ([7]) reads for 

r = r'^^ Rs{cos9,sm9), 9 = r^'', 

where we have used that i?s(r cos 6*, r sin 6*) = r''i?s (cos 6', sin 6^) since Rs{x,y) is a 
homogeneous polynomial of degree s. Hence, we see that the origin of system ([7]) 
is a focus if, and only if, the following integral is different from zero: 

/.27r 

£ = / i?, (cos^,sin^) d9 ^ 0. 
Jo 

We remark that if s is an odd number, then the origin of system ([7]) is a center. 
Easy computations show that 

Vo{x,y) = {x' + yY'^' 

is an inverse integrating factor for system ([7]) and that V{r, 9) = r*"'"^"^'^ is an 
inverse integrating factor of the ordinary differential equation on a cylinder corre- 
sponding to system ([7]). Thus, applying Theorem [U we deduce that the cyclicity 
Cjc[{Xi;,po) of the origin of system ([7]), with £ 7^ 0, is Cycl(A'£,po) > s/2. 

We observe that when A; = 0, we have that the origin of system ([7]) is a non- 
degenerate monodromic singular point and CycI{X^,pq) = s/2. The center problem 
for any even value of s is determined only by one Lyapunov constant, namely C, 
whereas the cyclicity of the origin (in case it is a focus) is given by s/2. We see, in 
this way, that the center problem and the determination of the cyclicity for a non- 
degenerate monodromic singular point are strongly related but are not equivalent 
problems. 

Example 5 In this example we show that no limit cycles bifurcate from a focus 
of a homogeneous system of degree d, under perturbations of subdegree > d. We 
consider a homogeneous system: 

X = Pd{x,y), y = qd{x,y), (8) 

where pd{x,y) and qd{x,y) are homogeneous polynomials of degree d with d > 1. 
A focus of system ([8]) has no characteristic directions because any linear real factor 
of ypd{x,y) — xqd{x,y) gives an invariant straight line of system ([8]) through the 
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origin. Since ypd{x, y) — xqd{x, y) has no linear real factors, we have that d is odd. 
In polar coordinates system ([H]) becomes 



e = f 



4-1 



where 



Pd{cos 9, sin 9) cos 9 + qd{cos 9, sin 9) sin 9, 
gd(cos 9, sin 9) cos 9 — pd(cos 9, sin 9) sin 9. 



The hypothesis that the origin is a focus implies that there are no characteristic 
directions. The non-existence of characteristic directions is equivalent to Fd{9) ^ 
for 9 G [0, 2tt). We see that the origin of system ([8]) is a focus if, and only if. 



On the other hand, easy computations show that Vo{x, y) = ypd{x, y)—xqd{x, y) 
is an inverse integrating factor for system ([8]), where we have used Euler's Theorem 
for homogeneous polynomials. Indeed V{r, 9) = r is an inverse integrating factor 
for the corresponding differential equation on the cylinder. By Theorem [1] we have 
that, when the origin of system ([H]) is a focus (i.e. £ 7^ 0), no limit cycles can 
bifurcate from it under perturbations of subdegree > d. 

However, if we take perturbations of lower degree, there can appear limit cycles 
which bifurcate from the origin. By Theorem [1], we have that at least {d — l)/2 
limit cycles bifurcate from the origin of system ([8]). 

For instance, let us consider the following system: 

X =ps{x,y) = {x - y){x'^ + y^) , y = q^{x,y) = {x + y){x'^ + y^), (9) 

which has an unstable degenerate focus at the origin. The function Vo(a;, 7/) = 
(x^ +2/^)^ is an inverse integrating factor for this system and, as we are under 
the hypothesis of Theorem [H we deduce that m = 1 and that no limit cycle can 
bifurcate from system ([9]) under perturbations of subdegree > 3. Let us consider 
the following perturbed system: 

x = (x- y){x^ + y'^) - e{x + y) , y = {x + y){x^ + y'^) + e{x - y), 

which has the invariant algebraic curve x"^ + y'^ = e. For e > 0, this algebraic curve 
is a hyperbolic limit cycle of the system which bifurcates from the origin. The 
function V^(x, = {x"^ + + y'^ — e) is an inverse integrating factor of this 

system. 




9 



1.2 A nilpotent focus 

We say that the origin of system ([1]) is a nilpotent singular point if it is a degenerate 
singularity that can be written as: 

X = y + P2{x,y) , y = Q2ix,y) , (10) 

with P2{x,y) and Q2{x,y) analytic functions near the origin without constant 
nor linear terms. The following theorem is due to Andreev [2] and it solves the 
monodromy problem for the origin of system ffTOl) . that is, it determines when the 
origin is a monodromic singular point. 

Theorem 6 [|2j| Let y = F{x) be the solution of y + P2{x,y) = passing through 
(0, 0). Define the functions f{x) = Q2{x, F{x)) = ax" + ■ ■ ■ with a ^ and a > 2 
and (j){x) = {dP2/dx + dQ2/ dy) {x , F{x)) . We have that either (pi^x) = hx^ + ■ ■ ■ 
with b ^ and f3 > 1 or (f){x) = 0. Then, the origin of (fTOj) is monodromic if, and 
only if, a < 0, a = 2n — 1 is an odd integer and one of the following conditions 
holds: 

(i) p>n-l. 

(ii) (3 = n — 1 and 6^ + Aan < 0. 

(iii) = 0. 

Definition 7 We consider a system of the form ffTOj) with the origin as a mon- 
odromic singular point. We define its Andreev number n > 2 as the corresponding 
integer value given in Theorem^ 

We consider system ( 1101) and we assume that the origin is a nilpotent mon- 
odromic singular point with Andreev number n. Then, the change of variables 

{x,y) (x,?/-F(x)), (11) 

where F[x) is defined in Theorem [HI and the scaling 

(x,y) ^ i^x,-^y), (12) 

with ^ = (— l/a)^/*^^~^"\ brings system (|T0|) into the following analytic form for 
monodromic nilpotent singularities 

X = y{-l+ Xi{x,y)), y = f{x) +y(f){x) +y'^Yo{x,y), (13) 

where Xi(0,0) = 0, f{x) = x'^"'~^ + ■■■ with n > 2 and either (f){x) = or 
4>{x) = bx^ + ■ ■ ■ with (3 > n—1. We remark that we have relabelled the functions 
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/(x), and the constant b with respect to the ones corresponding to system 
(fTUj) . We recall, cf. Theorem El that when /? = n — 1 we also have that b"^ — An < 0. 



We are going to transform system f|T3l) to an equation over a cylinder of the 
form ([3]). The transformation depends on the Andreev number n and it is given 
through the generalized trigonometric functions defined by Lyapunov [23j as the 
unique solution x{6) = Cs6 and y{6) = Sn^ of the following Cauchy problem 

^ = -y, I = x(0) = 1, y(0) = 0. (14) 

We observe that, in the particular case n = 1, the previous definition gives the 
classical trigonometric functions. 

We introduce in M^\{(0, 0)} the change to generalized polar coordinates, {x, y) ^— 
{r,9), defined by 

X = rCs9, y = r"Sn^. (15) 

In relation with this change, we say that a polynomial R{x, y) G C[x, y] is a (1, n)- 
quasihomogeneous polynomial of weighted degree w if the following identity is sat- 
isfied: 

RiXx, X-y) = rR{x,y), 

for all {x,y, A) G M^. We observe that a homogeneous polynomial of degree w is, 
with this definition, a (1, l)-quasihomogeneous polynomial of weighted degree w. 

Since we are going to use some properties of the generalized trigonometric 
functions and the relations satisfied by (1, n)-quasihomogeneous polynomials of 
weighted degree w, we summarize them up in the following proposition. The proof 
of each of its statements can be found in 



Proposition 8 [23| We fix an integer n > 1 and we consider {Cs6, Sn6) the 
solution of the Cauchy problem (HM . The following statements hold. 

Pk r (—) 

(a) The functions Cs9 and Sn9 are Tn-periodic with Tn = ^ y — (^"^+i^ ' 
where r(-) denotes the Euler Gamma function. 

(b) Cs^"0 + u'^vl'O = 1 (the fundamental relation). 

(c) Cs(-^) = CsO, Sn{-e) = -Sne, 

Cs{e + T.j2) = -Cse, Sn(0 + T„/2) = -Sn0. 

(d) Euler Theorem for quasihomogeneous polynomials: if R{x,y) is a (l,n)- 
quasihomogeneous polynomial of weighted degree w, then 

dR{x,y) dR{x,y) 
^ — dx — dy — ^ u!R{x,y)- 
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(e) Csifi = -Cse, Snip = (-l)"Sn^, where ip = (-1)"+^ + T„/2). // 



R{x,y) is a {1, n)-quasihomogeneous polynomial of weighted degree w, then 



R{-1,0) = i?(Cs(V2),Sn(V2)) = (-iri?(CsO,SnO) = (-iri?(l,0). 



Analogously to the case of a degenerate focus without characteristic directions, 
we can also provide the maximum number of limit cycles which can bifurcate from 
a nilpotent focus when only certain perturbations are taken into account. In this 
sense, and in relation with system ([2]), we consider the following definition, which 
will be used in the following Theorem [TOl 

Definition 9 An analytic perturbation vector field {P{x,y,6),Q{x,y,e)) is said 
to be (1, n)-quasihomogeneous of weighted subdegrees {wx,Wy) if P{\x,\^y,e) = 
C(A'^-) anc?_Q(Ax,A"y,£) = 0{X'"y). In this case, we denote by ;fi"'-"'"l = 
{P{x,y) + P{x,y,e))dx + {Q{x,y) + Q{x,y,e))dy the vector field associated to 
such a perturbation. 

We remark that the perturbative functions P{x,y,e),Q{x,y,e) do not need to 
be (1, n)-quasihomogeneous of a certain degree, they just need to have a (l,n)- 
quasihomogeneous sub degree high enough. 

The following theorem is one of the main results of this work. The symbol [xj 
denotes the integer part of x. 

Theorem 10 We assume that the origin of system (fTOl) is monodromic with An- 
dreev number n. Let V{r, 9) be an inverse integrating factor of the corresponding 
equation ^ which has a Laurent expansion in a neighborhood ofr = of the form 
V{r, 6) = v^{e) r"^ + 0(r'™+i), with Vm{0) ^0 andmeZ. 

(i) If m < or m + n is odd, then the origin of system fllOp is a center. 

(ii) // the origin of system (fTOj) is a focus, then m > 1, m + n is even and its 
cyclicity Cjcl^XsrjPo) satisfies Cycl(A'£,po) > (m + n)/2 — 1. In this case, 
m is the vanishing multiplicity ofV{r,6) onr = 0. 

(iii) // the origin of system (fT3l) is a focus and if only analytic perturbations of 
{1, n) -quasihomogeneous weighted subdegrees {w^^Wy) with > n and Wy > 
2n — 1 are taken into account, then the maximum number of limit cycles 



which bifurcate from the origin is [(m — l)/2j, that is, Cycl(A'£ " ,Po) = 



In particular. 




L(m-1)/2J. 
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The proof of this theorem is given in Section [21 



Remark 11 The proof of this theorem shows that if there exists an inverse in- 
tegrating factor Vq^XjH) of system (IT^ such that VQ{rCs9, Sn 9) / r'^"'^^ has a 
Laurent expansion in a neighborhood of r = 0, then the exponents of the lead- 
ing terms o/ Vo*(rCs6', r" Sn6')/r^"'~^ and V{r,9) coincide. Therefore, the value 
of m can he determined without performing the transformation of the system to 
generalized polar coordinates. 

The following corollary establishes a necessary condition for system ffTOl) to 
have an analytic inverse integrating factor Vo(x, defined in a neighborhood of 
the origin. 

Corollary 12 We assume that the origin of system ffTOj) is a nilpotent focus with 
Andreev number n, and that there exists an inverse integrating factor Vo{x,y) of 
flTU]) which is analytic in a neighborhood of the origin. Then, n is odd. 

The proofs of Theorem [10] and Corollary [12] are given in Section [3] Before the 
proofs of the main results of this paper we provide several examples of application 
of Theorem [TU] 

We would like to remark that the change to system ([T3|) is not always necessary 
to arrive at an equation over a cylinder, using generalized polar coordinates. The 
following proposition establishes sufficient conditions for an analytic system in 
cartesian coordinates (nilpotent or not) to be transformed to an equation over a 
cylinder by the change to generalized polar coordinates. We remark that given any 
analytic function P(x, y) in a neighborhood of the origin and any positive integer 
number n, we can always develop P{x, y) as a series of (1, ?7,)-quasihomogeneous 
polynomials. That is, we can always define (1, 'n.)-quasihomogeneous polynomials 
Pi{x, y) of weighted degree i such that the following identity is satisfied y) = 
Y.i>oVi{x,y). 

Proposition 13 Let h > 2 be an integer number and consider an analytic system 
of the form: 



where pi{x,y) and qi{x,y) are {l,n)-quasihomogeneous polynomials of weighted 
degree i and Pa{x, y), qb{x, y) are not identically null. Ifa — l = b — h>0 and 



for all 9 G [0,T^], then the origin of system (1161) is monodromic. Moreover, the 
change of coordinates x = rCs9, y = r"'Sn9 brings the system to an equation of 
the form ([3|) which is analytic in a neighborhood of its periodic solution r = 0. 





i>b 



(16) 



Cs9qb{Cs9,Sn9) -nSn9pa{Cs9,Sn9) ^ 
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Proof of Proposition [121 The change of variables {x, y) ^ (r, 9) gives 

. ^ x^^-^P{x,y)+yQ{x,y) _ r^^-^Cs^^"^g E.>a ^^P. W + ^^Sn g E.>6 ^^9^ W 

y,2fl—l ^2fj— 1 ' 

• ^ xQ{x,y)-fiyP{x,y) rCs9Y..^,r%{e) - fir^^neY^.^ypm 

where piiO) := pj(Cs6', Sn6'), gj(6') := gj(Cs6', Sn^). When a > 1 and b > h, this 
system is analytic at r = 0. In this case, we have that 

i>a i>b 

e = Cs^^r^-"g,(^) - nSn0^r^-^p,(0) . 

i>b i>a 

Moreover, we get that 

(Cs^ QbiO) - fi Sn9Pai9))r" 

where the dots denote terms of higher order in r. In this way, when a — 1 = b — n 
and Cs^ (lb{(^) — ^ Sn^ Pa{d) 7^ for all 9 G [0, Tj^), we have that 9 is different from 
zero for all |r| small enough. | 

We want to remark that not every system with a nilpotent singularity at the 
origin satisfies the hypothesis of Proposition[T3l Moreover, the value of the Andreev 
number n of a nilpotent system satisfying Proposition [13] does not need to coincide 
with the value of n which appears in the statement of this proposition. 





if 


a — 


1 > 6- 


h 




if 


a — 


1 < 6- 


h 


1 + - 


■ ■ if 


a — 


1 = 6- 


h 



Example 14 Let m and n be two positive integers such that n > 2, m > 1 and 
m + n is even. We consider the following system: 

X = y + xR{x,y), y = —x^^''^^ + ny R{x^y)^ (17) 

where y) is a (1, n)-quasihomogeneous polynomial of weighted degree m+n— 2. 
The origin of system (fTTll is a nilpotent singularity and it is a monodromic point 
applying Theorem [61 We observe that the Andreev number of system (ITTj) is n. 
The change to generalized polar coordinates (1151) brings the system to the form: 

f = r"+"~ii?(Cs0, Sn0), 9 = r"-\ 
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where we have used the fundamental relation Cs^^O + nSn^O = 1. Hence, the 
condition for the origin of system (fT7|l to be a focus is that 

C = R{Cse, SnO) de ^ 0, 
Jo 

where T„ is the period defined in Proposition [Hi For instance, if we take R{x, y) = 
always have that £ 7^ because the integrand is positive. The sym- 
metry properties of the generalized trigonometric functions Cs^^, SnO stated in 
Proposition [8] imply that in case m + n is odd, the value of C is zero. 
Easy computations show that 

m — 1 I 1 

Voix,y) = (x^" + ny')^'-' 

is an inverse integrating factor for system (ITTl) . where we have used Euler Theorem 
for (1, n)-quasihomogeneous polynomials, cf. Proposition [HI It is clear that since 
m and n are arbitrary positive integers with the restriction of being of the same 
parity, the function Vo{x, y) does not need to be analytic in a neighborhood of the 
origin. We observe that Vo{x, y) is analytic in a neighborhood of the origin if, and 
only if, there exists an integer k > such that m = 2kn + 1, which is an odd 
integer. We observe that: 

• When the origin of system (|T7|) is a focus, then m and n have the same parity. 

• When system ( 1T7I) has an analytic inverse integrating factor Vo{x,y) defined 
in a neighborhood of the origin, then m = 2kn + 1, which is an odd integer. 
Therefore, accordingly to Corollary [T2| we deduce that in this case and if the 
origin is a focus, then n needs to be an odd integer. 

The function V{r, 9) defined in the statement of Theorem [10] is V{r, 9) = 
which implies that the cyclicity Cycl(A'£,po) of the origin of system (ITTj) . when it 
is a focus, satisfies Cycl(A'e,po) > (m + n)/2 — 1. 

Another way to see that the origin of system ([TTIl is always a center in case 
that m + n is odd, is by using Theorem [TUl 

Example 15 We fix an integer n > 2 and we consider the following planar differ- 
ential system: 

X = y - uix"", y = -x^"-^ + z/2 x""^?/, (18) 

where i>i and z/2 are real parameters such that z/i z/2 — 1 < and (z/2 + n z/i)^ — 4n < 
0. We remark that the Andreev number of the nilpotent monodromic singular 
point at the origin of system f[T8l) is n. System f[T8|) is the most general (1, n)- 
quasihomogeneous planar polynomial differential system of weighted degree n with 
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a nilpotent monodromic singularity at the origin (we have used Andreev's Theorem 
E]). The change to generahzed polar coordinates (llSp transforms system (1181) in: 



e = -r"-i(l- (1^2 + 1^1^) Cs"eSn^). 

Hence, we have that the origin of system (fTS!) is a focus if, and only if, the following 
integral is different from zero: 

f^- Cs"-^g(-t/iCs^"g + t/2Sn^g) 
Jo 1 - (z/2 + 1^1 n) Cs^e Sue ^ ' 

where T„ is the minimal positive period of these generalized trigonometric func- 
tions, see Proposition [HI We define z{9) := 1 — (z/2 + n) Cs"6' Sn6' and we have 
that z{9) > for all ^ G M because 

z{e) = (^Cs"0 - ^l±^Sney + 1 (4n-(z/2 + riz/i)2)Sn2^, 

where we have used the fundamental relation Cs^^O + n Sn^^ = 1 and that, under 
our hypothesis, (z/2 + nvi)'^ — 4?t, < 0. 

We remark that, using the symmetries of the generalized trigonometric func- 
tions stated in Proposition [8], we deduce that if n is even, then C vanishes. There- 
fore, a necessary condition for the origin of system f|T8|) to be a focus is that n 
must be odd. Indeed, easy computations using the properties of the generalized 
trigonometric functions, show that 

Cs"-^^ (-Z/1 Cs^'^^ + z/2 Sn'^) _ z/2 - nz/i Cs""^^ ^ z\e) 



1- (z/2 + z/ir2)Cs"^Sne 2n z{e) 2n z{e) 

Hence, we deduce that 



2n Jq z{e) 

Under the assumptions n > 2 and odd, z/i z/2 — 1 < and (z/2 + nuiY — 4n < 0, 
we have that the integrand in the previous expression Cs^~^9/ z{9) > for any 
value of 6' G [0,T„]. Therefore, we have that, under these assumptions, the origin 
of system ( fTSl) is a focus if, and only if, z/2 ^ nui. 

The function Vo(x, y) = nyx — xy = x^" — (z/2 + z/i n) x^y + ny'^ is an inverse 
integrating factor of system ( JT8l) . which is a polynomial and, thus, it is analytic 
in the whole plane. We observe that the parity of n agrees with Corollary [T2l 
From the given expression of Vo(a;, y), we deduce that V{r,9) = r is an inverse 
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integrating factor of the equation corresponding to system (fTS]) . Applying 
Theorem [TUl we get that m = 1. Thus, if we consider an analytic perturbation 
which is (1, n)-quasihomogeneous of weighted subdegrees {wx,Wy), with Wx > n 
and Wy > 2n — 1, then no limit cycles can bifurcate from the origin. Indeed, we 
have that the cyclicity of the origin of system (flSj) is at least (n — l)/2. 

1.3 On the existence of an inverse integrating factor 

The following result is a summary and a generalization of several results on the 
existence of a smooth and non-flat inverse integrating factor Vo{x,y) in a neigh- 
borhood of an isolated singular point, see P| [T^ [2U] . 

Theorem 16 Let the origin be an isolated singular point of and let A,/i G C 
be the eigenvalues associated to the linear part of(^. If X ^ 0, then there exists 
a smooth and non-flat inverse integrating factor Vo{x,y) in a neighborhood of the 
origin. 

We can ensure the existence of an inverse integrating factor with stronger regu- 
larity in some particular cases. We recall that a singular point is said to be strong 
if the function dP/dx + dQ/dy is not zero on it and it is said to be weaA; otherwise. 

By using a translation, we can assume that the origin is the singular point under 
consideration. We say that the origin is analytically linearizable if there exists an 
analytic, near-identity change of variables such that the transformed vector field is 
linear. We say that the origin is orbitally analytically linearizable if there exists an 
analytic, near-identity change of variables such that the transformed vector field is 
a linear multiplied by a scalar unit function. 

Corollary 17 Let the origin be an isolated singular point o/dl]) and let A,/i G C 6e 
the eigenvalues associated to the linear part o/(II]). Then, the following statements 
hold. 

(i) (Strong focus) //A = a+iP and fi = a — i[3 with a, P G M\{0}, thenVQ{x,y) 
is analytic and it is unique up to a multiplicative constant. 

(ii) (Center) If X = i(3 and fi = —i(3 with (3 G IR\{0} and the origin is a center, 
then Vo(x, y) is analytic. 

(iii) (Linearizable point) // the origin is (orbitally) analytically linearizable and 
A 7^ 0, then Vo(a;, y) is analytic. 

(iv) (Node) // A,;U G M and A /i > 0, then Vo(a;,?/) is analytic and it is unique 
up to a multiplicative constant. 
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We include here the references where the previous statements have been proved. 
All the proofs are based upon the same idea: take the transformation to normal 
form, which is smooth in the considered cases, see |11J for instance. There is 
an analytic inverse integrating factor for the vector field in normal form (usually 
polynomial) and, thus, undoing the transformation, the obtained inverse integrat- 
ing factor is smooth. In this way, if the singular point is isolated and with A 7^ 0, 
then its (orbitally) linearizability implies the existence of an analytic inverse in- 
tegrating factor as it has been shown in [9]. If A 7^ 0, we have that the origin is 
either: 

- A strong focus, and the existence of an analytic inverse integrating factor is 
given in [201 El [13]. A strong focus is a particular case of a linearizable point. 

- A weak focus, that is, \ = i(3 and /i = —if3 with (3 G IR\{0} and the 
origin is not a center. The existence of a smooth inverse integrating factor is 
given in [T3], where a characterization of the existence of an analytic inverse 
integrating factor and an example of a weak focus without an analytic inverse 
integrating factor defined in a neighborhood of it are also given. In the 
following paragraph we include a sketch of the proof of this fact for the sake 
of completeness. 

- A non-degenerate center, and the existence of an analytic inverse integrat- 
ing factor is given by the normal form of a center given by Poincare. We 
recall that this normal form implies that a non-degenerate center is orbitally 
analytically linearizable. 

- A (hyperbolic) node, that is A,/i G M and A/i > 0, and the fact that the 
transformation to normal form is analytic is proved in [5j. Thus, there is an 
analytic inverse integrating factor defined in a neighborhood of it. Indeed, 
in a neighborhood of a node there cannot exist a first integral, see [9], so 
that the analytic inverse integrating factor is unique up to a multiplicative 
constant. 

- A (hyperbolic) saddle, that is A,/i G M and A/x < 0, and the existence of a 
smooth inverse integrating factor is shown in [T3] . 

- A semi-hyperbolic point, that is = but A 7^ 0. When it is isolated, the 
existence of a smooth inverse integrating factor is given in fO]. 

For the sake of completeness and in relation with the topic of this paper, we give 
the normal form in the case of a non-degenerate weak focus, that is, we consider 
the analytic system (|l]){f=o} and we suppose that its origin is a focus. In it 
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is shown the existence of a smooth and non-flat transformation that brings the 
considered system to the Birkhoff normal form: 



X = 



y + Si{r^)x - S2{r^)y 



y = X + S2{r'^)x + Si{r'^)y 



(19) 



where 5*1 and 5*2 are formal series on = + y"^. We use Borel's Theorem, see 
for instance [2Z] , to ensure the existence of smooth functions representing and 



Theorem 18 [Borel's Theorem] For every point p G M" and for every formal 
series in n variables, there exists a function f defined in a neighborhood of p 
whose Taylor series at p is equal to the given formal series. 

Therefore, we have a smooth change of coordinates which brings system (jl]){(=o} to 
system (1191) . Smooth changes of coordinates in normal form theory usually come 
from an order-to-order change and may induce flat terms in the normal form. 
However, since system (Il]){(=o} is analytic in a neighborhood of the origin, all of 
these flat terms can be removed by a suitable smooth change of coordinates. On 
the contrary, if system (ll]){^=o} were only smooth and with an infinite codimension 
focus at the origin (that is a focus with all its Liapunov constants equal to zero), 
then the normal form f|T9l) would be only formal because we cannot ensure the 
removal of all the flat terms. For example, the Birkhoff normal form of system 



is i; = — y = X but its origin is a focus. Therefore, the flat terms cannot be 
removed by any smooth change. This is an example of a smooth system with a 
focus of infinite codimension. It is known from Poincare that an analytic system 
(jlj) has no foci of infinite codimension. In this work we only consider analytic 
differential systems. 

We take the analytic system (|l]){(=o} and we perform the smooth change of 
coordinates which brings it to the Birkhoff normal form f|T9l) . It is well known, 
that the origin of system (|l]){^=o} is a center if, and only if. Si = 0. Since we are in 
the case that the origin is a focus, we have that 5*1 is not identically null and easy 
computations show that Vo{x,y) = (x^ + y^) S'i(x^ + y^) is an inverse integrating 
factor for the Birkhoff normal form. By undoing the change to the original system, 
we obtain a smooth and non-flat inverse integrating factor. 

For a strong focus, there is a unique analytic inverse integrating factor, as it has 
been shown in |9]. For a weak focus, when there is an analytic inverse integrating 
factor, then it is unique up to a multiplicative constant. However, there may exist 
other inverse integrating factors with lower regularity, as it is shown, for instance, 
in the forthcoming example with system fl22l) . 



^2. 
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The following result is a consequence of the results given in Theorems [T] and dU] 
and it ensures the existence of an inverse integrating factor, of class at least C^, in 
a neighborhood of certain degenerate focus points, namely for the origin of system 
([5]) without characteristic directions and the origin of system ffTOj) . 

Corollary 19 There exists an inverse integrating factor Vo{x,y), of class at least 
, in a neighborhood of the following two types of singular points: a degenerate 
focus without characteristic directions and a nilpotent focus. 

The proof of this corollary is given in Section [31 

2 Ordinary differential equations over a cylinder 



This section is devoted to several results related with ordinary differential equations 
of the form ([3]) defined over a cylinder C = {{r,6) E M. x : |r| < 5} for a certain 
S > sufficiently small. We denote by T the minimal positive period of the variable 
6, that is, we consider the circle = M/[0,T]. Thus, we consider an ordinary 
differential equation of the form Q: 

where J^{r, 9) is an analytic function on the cylinder C and 9) = 0. We have 
that, by assumption, the circle r = is a periodic orbit of equation Q. 

We assume that equation ([3]) has an inverse integrating factor V{r,9) which 
is analytic in a neighborhood of r = 0. All the results remain true if we assume 
that V{r, 9) is a smooth function in C which is non-fiat in r in a neighborhood of 
r = 0. We remark that V{r, 9) is a function over the cylinder C and, thus, V{r, 9) 
is T-periodic in 9. We recall that we have defined the vanishing multiplicity m of 
V{r, 6') on r = as the value such that 

V{r, 9) = Vm{9) r"" + O (r^+i) , (20) 

with Vm{9) ^ 0. The following lemma is already stated and proved in |15j. We 
include here a proof for the sake of completeness. 

Lemma 20 [15j IfV{r,9) is an inverse integrating factor of which is smooth 
and non-flat in r in a neighborhood of r = and with vanishing multiplicity m 
over r = 0, then the function Vm{9) defined in fl2Ul) satisfies that Vm{9) ^ for 
9 e [0,T). 
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Proof. By hypothesis J-'{0,6) = and we define the function Fi{6) as the one 
which satisfies 

J^{r,e) = Fi{e)r + O (r^) . 

We note that Fi{6) may be identically null. We have that V{r,6) satisfies the 
following partial differential equation 

and equating the coefficients of the order in the previous identity, we get that 
v'^{e) = (1 - m)Fi{e)vm{0). Since Vm{0) ^ 0, let G [0,r) be such that 
Vmifio) 7^ 0. We deduce that 



Vm{e) = vM) exp |(1 - m) F,{a) da j . (21) 
Therefore, we conclude that Vm{d) 7^ for G [0, T). | 



As we will see below, m coincides with the multiplicity of the limit cycle r = 
of equation ([3]). We remark that the integral -Fi(cr) da is the characteristic ex- 
ponent of the periodic orbit r = in equation ([3]). In particular, either m = 1 or 
Jq Fi{a) da = 0. Thus, from the formula (pTI) . we confirm that Vm{d) is always a 
T-periodic function. 

The following lemma establishes the existence and uniqueness of l^(r, 9) when 
the periodic orbit r = is isolated, that is, when r = is a limit cycle of ([3]). The 
existence is proved in [13] and the uniqueness is also stated and proved in [15] . 

Lemma 21 [13L 115] // the circle r = is a limit cycle of (I3l), then there exists 
an inverse integrating factor V{r, 6) of ([3]) which is smooth and non-fiat in r in 
a neighborhood o/r = 0. Indeed, V{r,9) is unique, up to a nonzero multiplicative 
constant. 

Proof. We recall here the main ideas to prove this statement, for the sake of 
completeness. 

Let us assume that r = is a limit cycle of multiplicity m of equation ([3]). 
When m = 1, we say that r = is hyperbolic. Following the ideas given in [13], 
and by a result in [30], in a neighborhood of r = 0, we can consider a smooth, 
and non-fiat in p in a neighborhood of p = 0, change of coordinates (r, 6) {p, r) 
which takes equation ([3]) to 

^ = X p with A 7^ 0, if m = 1; 

dr 

^ = + ap2'"~i with a G M, if m > 1. 
dr 
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In the case that r = is hyperbohc (m = 1), we have that the change of coordinates 
is, indeed, analytic in a neighborhood of r = 0. We remark that the function V{p, r) 
defined by: 

p if m = 1; 

p^ + ap^™-! if m>l; 

is an inverse integrating factor of the latter equation. By undoing the change of 
coordinates, we have a smooth, and non-flat in r in a neighborhood of r = 0, inverse 
integrating factor V^(r, 6) for equation ([3]). In the case that r = is hyperbolic, we 
have an analytic inverse integrating factor in a neighborhood of r = 0. 

To show the uniqueness of V{r,6), let us assume that there exist two linearly 
independent inverse integrating factors V{r, 9) and V{r, 9) of equation ([3]) which 
are both smooth and non-fiat in r in a neighborhood of r = 0. We assume that 
V{r,9) = v^{9)r"' + 0(r™+i) and V{r,9) = Vn^{9)r^ + 0(r'^+i) and that 
m > m. We have that the function on the cylinder C defined by H{r, 9) := 
V{r,9)/V{r,9) is not locally constant, smooth in r and of class in 9, by using 
the Lemma [201 

If m > m, we have that H{r, 9) is constant equal to all over the circle r = 
and if m = rh. using the proof of Lemma[20l we have that H(0. 9) = Vm{0) / Vm{0) , 
from which we deduce that H{r, 9) takes a constant value all over the circle r = 0. 
Moreover, this function H{r,9) is a first integral of equation ([3]), since it satisfies 
that 

o9 or 

Thus, H{r,9) is constant on each orbit of equation ([3]). If r = is a limit cycle, 
then the orbits in a neighborhood of this circle accumulate on it. By continuity of 
H{r, 9), this fact implies that H{r, 9) needs to take the same value on any point in 
a neighborhood of r = in contradiction with the fact that H{r, 9) is not locally 
constant. | 

The following example, which is described in page 219 of shows that the 
conditions for V{r, 9) to be smooth and non-fiat in r in a neighborhood of r = 
and 27r-periodic in 9 are essential to have a unique inverse integrating factor. We 
consider the planar differential system 

X = -y + x{x^ + y'^), y = x + yix"^ + y^), (22) 

which has a non-degenerate (unstable) focus at the origin. The following func- 
tions are two inverse integrating factors of the system which are of class in a 
neighborhood of the origin 

Vo{x, y) = (x^ + y'^Y and Vo{x, y) = (x^ + y'^Y sin ^2 arctan ^— j + 



x^ + y^ 
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In polar coordinates, this system reads for r = r^, 9 = 1, which has the follow- 
ing inverse integrating factor Vi{r,6) = Vo(r cos^,rsin^)/r = analytic in a 
neighborhood of r = 0. Moreover, the function V2{r,6) = Vo{r cos6,r sm6)/r = 
sin(2^ + r~^) is an inverse integrating factor of class in r 7^ 0. Indeed, 
V3{r,6) = r + 2r^6 is another inverse integrating factor of the system in polar 
coordinates, which is analytic in r but not 27r-periodic in 6. 

We will show that when the periodic orbit r = is a limit cycle of equation 
([3]), the vanishing multiplicity of an inverse integrating factor V{r,6) is strictly 
positive. 

Lemma 22 Let us consider a differential equation of the form ([3]) over a cylinder 
C = {{r,6) E X : |r| < 5} for a certain 6 > and let us assume that r = 
is a periodic orbit of the equation. We assume that V{r, 9) is an inverse integrating 
factor for equation ([3]) defined in C\{r = 0} and which has a Laurent series of the 
form: 

V{r,e) = v^{e)r"' + C(r'"+i) , 

with Vm{0) ^ and m G Z. If m < 0, the periodic orbit r = has a neighborhood 
filled with periodic orbits, that is, it is not a limit cycle. 

Proof. Let u) = dr — JF(r, 6) dO be the Pfaffian 1-form associated to equation 
([3]). Since V{r.,d) is an inverse integrating factor of the equation ([3]), we have 
that uj /V is a closed 1-form. We observe that the proof of Lemma [2U] also applies 
and, therefore, we have that Vm{6) 7^ for any 6 G [0,T). In case that m < 0, 
we have that u/V is well-defined in the whole cylinder C. Let us consider any 
non-contractible cycle in this cylinder, for instance the cycle r = 0. By virtue 
of De Rham's Theorem, see [H], we have that the closed 1-form u/V is exact 
if, and only if, the value of the line integral J^^qUj/V is zero. This value is zero 
since the oval r = is an orbit of the 1-form u and, thus, co'|r-=o = 0. Hence, 
we have that u/V = dH for a certain function H{r,6), which turns out to 
be a first integral of equation ([3]). The existence of this first integral implies that 
the cycle r = is surrounded by periodic orbits, formed by the level curves of H. | 

In order to relate the vanishing multiplicity m of V{r, 6) on r = 0, which we 
have proved to be positive, and the cychcity of the focus at the origin of system 
([1]), we use a previous result which is already stated and proved in [15]. Our result 
gives an ordinary differential equation for the Poincare map associated to equation 
([3]) in terms of the inverse integrating factor V{r, 9). If the minimal positive period 
of ^ in ([2]) is denoted by T and \E'(^;ro) is the flow of ([3]) with initial condition 
^(Os'^o) = ''"o; "we recall that the Poincare map 11 : S C M ^ M is defined as 
n(ro) = ^'(T;ro). We have that 11 is an analytic diffeomorphism defined in a 
neighborhood S of tq = whose fixed points correspond to periodic orbits of the 
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equation ([3]). Since we have that r = is a hmit cycle of equation ([2]), we deduce 
that the Poincare map is not the identity map. We recall that the limit cycle r = 
is said to have multiplicity k if the expansion of the analytic Poincare map in a 
neighborhood of tq = is of the form n(ro) = tq + Tq + 0{r^'^^), where Ck ^ 0. 
The multiplicity of the limit cycle r = in equation ([3]) allows us to determine the 
cyclicity of the focus at the origin of system ([T]), as we will see below. 

We state the result proved in [12] in the terms used in the present paper. 



Theorem 23 [15| Let us assume that V{r, 9) is an inverse integrating factor of 



equation ^ defined in a neighborhood of the periodic orbit r = 0, whose minimal 
positive period is denoted by T. We consider IKjq) the Poincare map associated to 
the periodic orbit r = of equation ([3]) . Then, the following identity holds: 

V{Il{ro),T) = V{ro,0) n'(ro). (23) 

As a consequence of equation (!23|) . we can prove the following result. 

Corollary 24 [15] Let us assume that V{r,6) is an inverse integrating factor of 
equation ([3]) which is smooth and non-flat in r in a neighborhood of the limit cycle 
r = and whose vanishing multiplicity over it is m. Then, r = is a limit cycle 
of multiplicity m. 

Proof. Since V{r, 9) is assumed to be a function on the cylinder C, we have that it 
needs to be T periodic in 0. We have that V{rQ, T) = V{rQ, 0) and we consider its 
development in a neighborhood of Tq = 0: V{rQ,0) = Umr^ + O {r^~^^), where 
I'm 0. We observe that the index m appearing in this decomposition coincides 
with the vanishing multiplicity of V{r,6) in r = by Lemma [201 Recalling that 
n(ro) = To + CfcTQ + 0{rQ~^^), where Ck ^ 0, we consider equation fl25]) and we 
subtract V {tq, 0) from both members to obtain that: 

The lowest order terms in both sides of the previous identity correspond to r^^^~^ 
and the equation of their coefficients is mcki'm = kckVm-, which implies that 
k = m. I 



3 Proofs of the results 
Proof of Theorem [B 

The following lemma establishes the first step in the proof of Theorem [TJ We show 
that the transformation to polar coordinates x = rcos6, y = rsin6', of system ([5]) 
gives an equation over a cylinder of the form ([3]). 
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Lemma 25 We consider system with d > 1 and odd integer and we assume 
there are no characteristic directions. Then the transformation to polar coordinates 
brings system ([5]) to an ordinary differential equation over a cylinder. 



(24) 



Proof. In polar coordinates system ([5]) becomes 

r = r'^R{r,9) = r'^ {Rd{9) + O (r)) , 
Q = r^"iF(r,^) = r'^-i(Fd(^) + 0(r)), 

where 

-Rd(^) = Pd(cos^^, sin^^) cos6' + grf(cos6', sin6') sin6', 
Pd{&) = gd(cos6', sin 6') cos 6* — prf(cos6', sin6') sin6'. 

The hypothesis that there are no characteristic directions is equivalent to say that 
7^ fo^^ ^ ^ [0, 27r). We can, therefore, consider the ordinary differential 
equation associated to the orbits of system fl2^ which takes the form ([3]): 

^ ^ rR{r,e) 

de F{r,e) ■ ^ ' 

■ 

The following lemma establishes that the center problem for the origin of system 
([5]) is equivalent to the determine when the circle r = is contained in a period 
annulus for equation ( !25l) . 

Lemma 26 We consider system ([5]) with d > 1 and odd integer and we assume 
there are no characteristic directions. The circle r = is a limit cycle of equation 
(!25|) if, and only if, the origin of system is a focus. 

Proof. The origin of system ([5]) is transformed to the periodic orbit r = of equa- 
tion (125|) by the transformation to polar coordinates. This transformation gives a 
one-to-one correspondence between any point in a punctured neighborhood of the 
origin in the plane {x, y) and a cylinder {{r,6) : < r < 6,6 & S^} for 6 > suf- 
ficiently small. Thus, any orbit spiraling from or towards the origin of system ([5]) 
is transformed to an orbit spiraling (from or towards) the circle r = in equation 

m- ■ 

We have a symmetry for equation (l25l) which is inherited by the symmetries of 
the polar coordinates. 

Lemma 27 Let us consider a planar differential system x = P{x,y), y = 
Q{x,y), and perform the change to polar coordinates x = r cos6, y = r sinO. 
The resulting system r = 'E{r,9), 9 = Q{r,6) is invariant under the change of 
variables (r, 6) i— > (— r, 6 + tt) . 
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Proof. We observe that the monomials r cos 6 and r sin 6 are invariant by the 
change of variables. Since 



e = 

the result follows. 



r cos 6 P{r cos 6*, r sin 6*) + r sin 9Q{r cos ^, r sin i 

r 

r cos 9Q{r cos 0, r sin 9) ~ r sin ^ P(r cos 9, r sin i 



Lemma 28 We assume that the origin of system is a focus without charac- 
teristic direction and we consider the corresponding equation ( !25|) with associated 
Poincare map n(ro) = tq + c^r'^ + with c„ ^ 0. Then, m is odd. 

Proof. By Lemma [26], we have that the origin of system ([5]) is a focus if, and only 
if, the circle r = is a limit cycle of equation fl25l) . We assume that the circle 
r = is a limit cycle with multiplicity m and we consider its associated Poincare 
map n(ro). By Lemma [271 we have that equation f[25|) has the discrete symmetry 
(r, ^) hH> {—r,9 + tt) because it comes from a system in cartesian coordinates 
This symmetry implies that r = is either a stable or an unstable limit cycle. 
Thus, m needs to be odd. | 

The next lemma states that, under our hypothesis, we have an inverse inte- 
grating factor V{r, 9) for equation (|25l) which is smooth and non-fiat in r in a 
neighborhood of r = and such that V^(0, 9) = 0. 

Lemma 29 We consider system ([5]) with d > 1 and odd integer and we assume 
there are no characteristic directions. 

(i) If system ([5]) has an inverse integrating factor Vo{x , y) defined in a neighbor- 
hood of the origin, then the function defined by 

Vo{rcos9,rsm9) 

V[r,9) := 



rdF{r,9) 

is an inverse integrating factor for equation ( l25i) in r ^ 0. 

(ii) Let V{r, 9) be an inverse integrating factor of equation f[25l) which has a Lau- 
rent expansion in a neighborhood o/r = of the form V{r, 9) = Vm{9) r"^ + 
0{r"^~^^), with Vm{9) ^ and m G Z. Then, if m < the origin of system 
([5]) zs a center. 
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Proof, (i) Since the jacobian to polar coordinates is r, we have that the func- 
tion Vo(rcos^, r sin^)/r is an inverse integrating factor for system (1241) . We see 
that equation fl25|) is the equation of the orbits associated to system fl2^ and, 
therefore, the function V{r,6) := Vo(r cos 6', r sin 6')/(r'^ F(r, ^)) is an inverse inte- 
grating factor of fl25|) . We remark that V{r, 6) does not need to be well-defined in 
a neighborhood of r = 0. Thus, we only have, at the moment, that it is an inverse 
integrating factor for equation fl25l) in r 7^ 0. 

(ii) If m < 0, we are under the hypothesis of Lemma [221 and we conclude that 
the cycle r = is surrounded by periodic orbits, which give rise to a neighborhood 
of the origin of system ([5]) filled with periodic orbits. Therefore, the origin of sys- 
tem ([S]) is a center. | 

The previous Lemma [29l together with the statements given in Lemmas [25] 
and [26], establishes that under the hypothesis of Theorem [T] we have that m > 1. 
Moreover, by Lemma [28] and Corollary [2ll we have that m needs to be odd. We 
have proved statement (i) in Theorem [1] 

Assuming that the origin of system ([S]) is a focus, the following step of the 
proof is to relate the multiplicity of the limit cycle r = of equation (!25|) and the 
cyclicity of the origin of system ([5]). 

Lemma 30 We consider system (jSj) with d>l and odd integer, we assume there 
are no characteristic directions and that the origin po is a focus with cyclicity 
Cycl(A'£,po)- We consider the corresponding equation (1251) and we assume that 
r = is a limit cycle of multiplicity m = 2/c+l. Then, Cjcl{Xe,po) > {m+d)/2 — l. 
Moreover: 

1. When d = 1, we have Cycl(A'£,po) = k. 

2. // we only consider perturbations of system whose subdegree is > d, then 
at most k limit cycles can bifurcate from the origin of system that is, 

Cycl(A^i"],Po) = k. 

Proof. The particular case d = 1 of statement 1. and system ([5]) with a focus at 
the origin is proved in Theorem 40, Ch. IX in [4] (page 254), see also [29] . 

We first provide an example of a perturbation of equation ([25]) . with m limit 
cycles bifurcating from r = 0, whose transformation to cartesian coordinates gives 
a perturbation of system (jSj) of the form with exactly (m — 1)/2 = k limit cycles 
bifurcating from the origin. This example shows that Cyc\{X^,pq) > k. 

We take the corresponding equation f[25l) and we assume that r = is a limit 
cycle of multiplicity m, which is an odd integer with m > 1. We consider the 
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associated system fl24l) from which fl25p comes from, we define A; = (m — l)/2 and 
we perturb system fl24p in the following way: 

k-l 

r = r'^ R{r,e) + ^e^-' air'^'+\ 9 = r'^-^ F{r,e), (26) 

1=0 

with the convention that if = no perturbation term is taken. The real constant 
e is the perturbation parameter < \e\ « 1 and the a^, = 0, 1, 2, . . . , — 1, are 
real constants to be chosen in such a way that the Poincare map 11 (ro; e) associated 
to the ordinary differential equation 

dr _ r'^ R{r, 9) + J2i=i _ ^^(^, 9) + r Y!IZq e'"'" ca r'^' 

d9 ~ ^d-i F(^r, 9) ~ F{r,9) 

has 2k + 1 real zeroes; k of them positive. We recall that F{0, 9) ^ for all 
9 G [0, 2tt) and, thus, the perturbative terms are an analytic perturbation in a 
neighborhood of r = and e = of equation fl25l) . The proof of the fact that this 
choice of Oj can be done is analogous to the one described in [4j, pp 254-259. More 
concretely, the exponent of the leading term of the displacement function d{rQ; 0) 
of system fl25l) is m and the considered perturbation fl26l) produces that d{ro; e) 
has all the monomials of odd powers of tq up to order m. The coefficient of each 
monomial, for e sufficiently small, is dominated by one of the constants Oj. 

Undoing the change to polar coordinates, system fl2B]) gives rise to an analytic 
system in a neighborhood of the origin which is a perturbation of system of 
the form ([2]) and with A; = (m — l)/2 limit cycles bifurcating from the origin. If 
system (|5]) is written as a; = P{x, y) and y = Q{x, y), then the change to cartesian 
coordinates from (l26l) reads for: 

X = P{x,y) + xK{x,y,e), y = Q{x,y) + yK{x,y,e), (27) 
fc-i 

where K{x,y,€) = ^£:'^~*ai (x^ + y'^y^~ . We recall that d is odd and d > 1. 

2 = 

In this way, we have that Cyd{X^,po) > k. 

We provide now an example of an analytic perturbation of system ([5]) with 
at least {m + d)/2 — 1 limit cycles bifurcating from the origin. We take system 
f l27|) and we perturb it in order to produce i = (d — l)/2 additional limit cycles 
bifurcating from the origin when e — 0. Let us consider the smallest limit cycle 
7 surrounding the origin of system (p7|) and let us assume that it is an attractor. 
We have that the (m — l)/2 limit cycles of system (127|) which bifurcate from the 
origin are hyperbolic, by choosing the parameters Oj conveniently. Since 7 is the 
smallest limit cycle, we have that the origin needs to be a repeller. Let us take a 
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convenient real value such that the system 

X = P{x,y) + X K{x,y,e) + e'''^^ xbe_i{x'^ + y'^Y'^, 
y = Q{x, y) + y K{x, y, e) + e^'^^ y (x^ + y^-' 

still has the previous limit cycle 7 as an attractor (thus, \bg^i\ needs to be small 
enough) and the origin also becomes an attractor (this implies that bi^i < 0). 
Therefore, there is a limit cycle bifurcating from the origin, surrounding it and 
smaller than 7. If 7 was a repeller and the origin an attractor in system fl27j) . we 
take bi^i > in order to make the new limit cycle to bifurcate. This bifurcated limit 
cycle is hyperbolic by conveniently choosing the value The previous system 
maintains the (m — l)/2 limit cycles of flTTI) because they are all hyperbolic. Thus, 
the previous system has at least (m — l)/2 + 1 limit cycles bifurcating from the 
origin when e — >■ 0. 

By induction, and a relabelling of the parameters at := fei+(d-i)/2, "we deduce 
that the following system 

X = P{x,y) + xK{x,y,e), y = Q{x,y) + yK{x,y,e), (28) 

L-l 

where K{x,y,e) = ^ e^^"* 6j (x^ + y'^f, L := (m + d)/2 - 1, has at least 

{m + d)/2 — 1 limit cycles bifurcating from the origin. We recall that both m and 
d are odd and d > 1, m > 1. In this way, we have that Cycl(A'e,po) ^ (m + d) /2 — 1. 

Finally, we will prove statement 2. If we only consider perturbations of system 
dSD whose subdegree is > d, that is vector fields of the type A'J'^' , then, the transfor- 
mation of these perturbative terms to polar coordinates gives rise to a perturbation 
of the corresponding equation fl25|) which is analytic in a neighborhood of r = 
and e = 0. Let us assume that the circle r = is a limit cycle with multiplicity 
m and we consider the Poincare map n(ro) defined in the previous section, which 
satisfies: 

n(ro) = ro + c^C + 0{r^^'), 

with Cm 7^ 0. We recall that equation fl25]) has the discrete symmetry (r, 6) 
{—r,6 + tt) because it comes from a system in cartesian coordinates ([5]). This 
symmetry implies that r = is a limit cycle which cannot be semistable and 
therefore m is odd. 

The key point of the proof is that any perturbation of (!25|) analytic near (r, e) = 
(0, 0), with e G small, has a displacement function d{rQ; e) = n(ro; e) — tq which 
is analytic near (ro, e) = (0, 0) and when e = coincides with the displacement 
function d{rQ] 0) = n(ro) — ro of the unperturbed equation fl25l) . By using standard 
arguments (counting zeroes with Weierstrass Preparation Theorem of d{ro] e) near 
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{rQ,e) = (0,0)), the cyclicity of the circle r = under analytic perturbations of 
equation (1251) is m. We recall that the cyclicity and the multiplicity of a limit cycle 
are equal, see [3]. However, since the displacement function d{ro; 0) = n(ro) — ro 
of the unperturbed equation fl25l) is of odd order m at ro = 0, and taking into 
account the above mentioned discrete symmetry, we have that only (m — l)/2 
zeroes of d{rQ] e) can appear for tq > and ||£:|| small enough. This fact gives that 
at most (m — l)/2 limit cycles bifurcate from the origin pq of system when only 
this kind of perturbative terms are taken into account. Therefore, we have proved 
that CycI{X£,po) = (m — l)/2. Indeed, the example given in (1271) shows that this 
upper bound is sharp. | 



Proof of Theorem [TUl 

The proof of this theorem goes analogously to the proof of Theorem [H only with 
some technical differences. 



We first show that the transformation to generalized polar coordinates fllSp 
X = rCs6,y = r^SnO transforms system f|T3l) to an equation over the cylinder 
of the form ([3]). 

Lemma 31 We assume that the origin of system (1131) is a nilpotent monodromic 
singular point. Then the transformation to generalized polar coordinates {x,y) = 
{rCs9, r"'Sia9) brings system ( IT3l) to an ordinary differential equation ([3]) over a 
cylinder. 

Proof. Taking into account that Cs^^O + nSn^9 = 1, we get that the Jacobian 
determinant of the former change is 

' ^ d{r,e) 
Since x^"" + ny"^ = r^", we deduce that 

^2n-i yy xy — nyx 

T = u = . 

j,2n—l ' joji+l 

In particular, system (fT3!) adopts the form 

r = p(^)r"+i + 0(r"+2)^ = ^'^-i +0(r'^), (29) 

when j3 > n — 1 or (j){x) = 0, and 

f = 6Cs"-^^Sn20r" + O(r'^+^), ^ = (1 + 6 Cs"^ Sn 0) r"-^ + 0(r"), (30) 
when (3 = n — 1. 
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We also observe that in the latter case (/? = n — 1) we have the following 
decomposition 



1 + b Cs"^ Sn ^ = (^Cs'^e + ^SnO^ +^{An- 6^) Sn^^, 

where we have used that Cs^"^ + n Sn^O = 1. Since 4?t, — 6^ > in this case, 
due to Andreev's conditions for monodromy, we deduce that 1+6 Cs"^ Sn 6' > 
for any ^ G M. 

We denote by E{r,9) the function defined by r and by <d{r,9) the function 
defined by 6 in both cases, and we have that, 

r = E{r,e), 9 = Q{r,e) = e„_i(^)r"-i + O (r") , (31) 

where 

1 if /? > n - 1 or = 0, 

1 + 6Cs"^Sn^ if (3 = n-l. 

Hence, the equation of the orbits corresponding to system (l29l) or (130|1 writes as 

0(r2) 



dr 



1 + 0(r) 

0{r) 



l + 6Cs"^Sn0 + O(r) 



if (3 > n — 1 or (f){x) = 0, 
if (3 = n-l. 



We observe that 9„_i(6') ^ for any 9 G [0, T„). In short, we have proved that in a 
neighborhood of any monodromic singular point of the form ffTOl) . we can perform 
a transformation, which is the composition of the changes (|TT1) and (|T2l) and the 
transformation to generalized polar coordinates, which brings the system to an 
equation over a cylinder C: 

To = ^"••«)- 

where jF(r, 9) is T„-periodic in 9 and J-'{0, 9) = 0. | 

The center problem for the origin of system (ITOl) is equivalent to determine 
when the circle r = is contained in a period annulus for equation (|32l) . 



Lemma 32 VKe assume that the origin of system (11 01) zs a monodromic singular 
point with Andreev number n. The circle r = is a limit cycle of equation fl32|) if, 
and only if, the origin of system ( fTOl) is a focus. 
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Proof. The proof is analogous to the proof of Lemma [2B The transformation from 
system ( TTUj) to equation (1521) gives a one-to-one correspondence between each point 
in a punctured neighborhood of the origin in the plane {x, y) and each point on a 
cylinder {{r,6) : < r < 6,6 E S^} for 6 > sufficiently small. | 

The following proposition establishes a symmetry for equation 0321] which is 
inherited by the symmetries of the generalized trigonometric functions. 

Proposition 33 Let us consider a planar differential system x = P{x,y), 
y = Q{x,y), we take any positive integer n and perform the change to generalized 
polar coordinates x = rCs6, y = r"Sn^. The resulting system r = 'E{r,9), 9 = 
Q{r, 6) is invariant under the change of variables [r, 6) i— > (— r, (— l)""*"^ [6 + Tn/2]). 

Proof. We observe that, due to Proposition [HI the following composition is the 
identity {X = x, Y = y): 

ix,y) ^ {r,9) ^ (R,^) ^ {X,Y), 

where R=-r,ip= + r„/2] and X = RCsif, Y = Sinp. Since 

. ^ x^''-^P{x,y) + yQjx^y) ■ ^ xQ{x,y) - nyP{x,y) 

the proposition follows. | 

The previous symmetry of equation (J32i) imposes a condition on the circle r = 
to be a limit cycle. 

Lemma 34 We assume that the origin of system fllOl) is a focus with Andreev 
number n. We consider the corresponding equation ( 1521) and its Poincare map 
n(ro) = ro + c„C + 0{r^^^), with Cm ^ 0. 

(i) If n is odd, then r = cannot be a semistable limit cycle of equation (!32|) . 
that is, m is odd. 

(ii) If n is even, then r = is a semistable limit cycle of equation (!32l) . that is, 
m is even. 

Proof. The change to generalized polar coordinates ensures that the origin of 
system (fTOj) is a focus if, and only if, the circle r = is a limit cycle of equation 
(!32|) . We assume that the circle r = is a limit cycle with multiplicity m and 
we consider its associated Poincare map n(ro). We recall that equation (152]) has 
the discrete symmetry {r,6) i— > (— r, (—1)"+^ [6 + r„/2]) because it comes from a 
system in cartesian coordinates (|TOl) . see Proposition [331 This symmetry implies 
that r = is a semistable limit cycle if, and only if, n is even. | 
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Corollary 35 // equation 0321) has a periodic orbit different from r = 0, then it 
has two periodic orbits (one in the upper half cylinder and one in the lower half 
cylinder). 

Proof. The properties of the periodic orbits of fl32|) stated in this corollary are 
straightforward consequences of the discrete symmetry given in Proposition [331 | 

Let Vo{x,y) be an inverse integrating factor defined in a neighborhood of the 
origin for system (fTOl) . Since the changes of variables (fTT!) and (fT2|) have con- 
stant non-vanishing Jacobian, it follows that the transformed system ( [T3|) has the 
following inverse integrating factor defined in a neighborhood of the origin: 

V*ix,y) = Voir'x,-r'y + Fir'x)). 

Next lemma is the analogous to Lemma [291 and states that, under our hypothe- 
sis, we have an inverse integrating factor V{r, 9) for equation (1321) which is analytic 
in r in a neighborhood of r = and such that V^(0, 9) = 0. 



Lemma 36 We assume that the origin of system (1131) is a nilpotent monodromic 
singularity. 

(i) If system (IT^ has an inverse integrating factor Vq*(x, ?/) defined in a neigh- 
borhood of the origin, then the function defined by 

V^r^e) := VoirCs9,r-Sn9) 



r"e(r,^) 

where 0(r, 9) is the function defined in (^T\i . is an inverse integrating factor 
for equation ( H^ in r 0. 

(ii) Let V{r,9) be an inverse integrating factor of equation f[32l) . We assume 
that V{r, 9) has a Laurent expansion in a neighborhood of r = of the form 
V{r, 9) = v^{9) + C»(r"*+^), with Vmi9) ^ and m e Z. Then, zf m < 
the origin of system (lT3l) is a center. 

Proof, (i) Taking into account the Jacobian of the change to generalized polar 
coordinates {x,y) = {rCs9, r^Sn9), we see that the differential equation (l32l) 
has the inverse integrating factor V{r, 9) described in the statement which is a 
T„-periodic function of 9. We observe that V{r, 9) may not be well-defined on 
r = 0. 

(ii) The assumption m < establishes that we are under the hypothesis of 
Lemma [221 and we conclude that the cycle r = is surrounded by periodic orbits, 
which give rise to a neighborhood of the origin of system ( |T3l) filled with periodic 
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orbits. 
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The previous Lemmas [Ml and [36| together with Corollary [241 ensure that, un- 
der the hypothesis of Theorem [TOl and if the origin of system (ITOl) is a focus, then 
m > 1 and m + n needs to be even. Thus, we have proved statement (i) in Theorem 

m 

To end with, we relate the multiplicity of the limit cycle r = of equation (!32|) 
and the cyclicity of the origin of system ([TOj) . 

Lemma 37 We assume that the origin of system f [TO|) is a nilpotent focus with 
Andreev number n. We consider the corresponding equation ( !32l) for which we 
assume the circle r = to be a limit cycle with multiplicity m. 

1. The cyclicity Cjcl{Xs,po) of the origin of system ( JTOjl satisfies Cjc\{Xs,po) > 
(m + n)/2- 1. 

2. If only analytic perturbations of system (113!) with {l,n)-quasihomogeneous 
weighted subdegrees {wx,Wy) with > n and Wy > 2n — 1 are taken into 
account, then the maximum number of limit cycles which bifurcate from the 
origin of system ([T51) {and, thus, of system ( TTUl) ) is [(m — l)/2j, that is, 
CYd{xt'''-'\po) = [{m-l)/2\. 

Proof. The fact that m and n have the same parity is a consequence of Lemma [Ml 
If n is odd, the symmetry given in Proposition l33l implies that at most (m — l)/2 
limit cycles can bifurcate from the limit cycle r = of equation f[32|) in the region 
r > 0. If is even, since equation f[32l) comes from cartesian coordinates, we also 
need to take into account that r = is always a solution. Therefore, by using the 
symmetry again, at most (m — 2)/2 limit cycles can bifurcate from r = in the 
region r > 0. 

We provide an example of a perturbation of equation (1321) . with m limit cycles 
bifurcating from r = (counting multiplicities), whose transformation to cartesian 
coordinates gives a perturbation of system (fTOj) with exactly [(m — 1)/2J limit 
cycles bifurcating from the origin. This example proves that Cyc\{X^,pq) > k. We 
take equation (1321) and we assume that r = is a limit cycle of multiplicity m, 
which is an integer with the same parity as n and such that m > 1. We consider 
the associated system f[3Tl) from which f[32|) comes from, we define k = \_{m — 1)/2J 
and we perturb system f[3T|) in the following way: 

k-1 

r = S(r, 9) + Y, (Cs ^ g ^ q^^^ ^ ^^^^ 

T-i (33) 

f = S(r,e) + ^e'=-^air"+^+2^(Cs5)"+^% ^ = e(r,0), if n is even, 

j=0 
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with the convention that if A; = no perturbation terms are taken. The real value 
e is the perturbation parameter < |£:| << 1 and the a^, i = 0, 1, 2, . . . , — 1, 
are real constants to be chosen so that the Poincare map n(ro; e) associated to the 
ordinary differential equation 



dr 



a,- r 



n+2i 



(Cs( 



\n-l+2i 



e(r,^) 



„n+l+2i 



(Cs^) 



n+2i 



e{r,e) 



if n is odd, 
if n is even, 



has 2k + 1 real zeroes; k of them positive. The proof of this fact is analogous to 
the proof given in pp. 254-259. 

Undoing the change to generalized polar coordinates system f l33|) gives rise to 
an analytic system in a neighborhood of the origin which is a perturbed system 
from f|T3|) of the form ([2]) and with k = \_{m — l)/2j limit cycles bifurcating from 
the origin. All these limit cycles are hyperbolic by taking convenient values of the 
parameters a^. If system (IT^ is written as i = P{x,y) and y = Q{x,y), then the 
change to cartesian coordinates from (133!) reads for: 



X 



P{x,y) + xK{x,y,e), y = Q{x,y) +nyK{x,y,e) 



(34) 



where 



K{x,y,e) 



fe-i 

E 

i=0 



n+2i 



if n is odd, 



if n is even. 



Then, undoing the changes ffTTl) and f|T2l) . we obtain a perturbation of system ffTOj) 
with = [(m — l)/2j limit cycles bifurcating from the origin. We observe that 
these perturbative terms satisfy the (1, ?7,)-quasihomogeneous subdegree conditions 
established in the statement 2 of the lemma. 

In order to generate the i := \n/2\ limit cycles that we lack to prove the first 
statement of this lemma, we will consider perturbations with lower subdegree. We 
observe that the stability of the origin in system fl34l) is given by the sign of Oq. 
Thus, if we choose a convenient value be-i such that the following system 



.2{<?-l) 



X = P{x,y) + X K{x,y,e) + e xbe^ix 
y = Q{x,y) +nyK{x,y,e) + e'^'^^ybi^ix'^'^^'^^ 

satisfies that the stability of the smallest limit cycle in fl34|) does not change (thus 
needs to be small enough) and with the origin of the contrary stability 
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(namely, of contrary sign to Oq), then there is a new hmit cycle bifurcat- 
ing from the origin. This limit cycle is hyperbolic by choosing 6^„i conveniently. 
By induction, and relabelling := we have that the system 

X = P{x,y) + xK{x,y,e), y = Q{x,y) +nyK{x,y,e), (35) 

where 

L-l 

Kix,y,e) = 

i=0 

and L = (m + n)/2 — 1, has at least (m + n)/2 — 1 limit cycles bifurcating from 
the origin. By undoing the changes (fTT]) and (fT2l) . we have thus proved that the 
cyclicity of the origin of system (fTOj) is at least {m + n)/2 — 1. 

We are going to prove statement 2 of this lemma. By the results established 
on Section [2l we can control the maximum number of limit cycles which bifurcate 
from r = in equation fl32|) under analytic perturbations of this equation. The 
hypothesis that we only take analytic perturbations of system f fT3|) with {l,n)- 
quasihomogeneous weighted subdegrees (wxyWy) with > n and Wy > 2n — 1 is 
equivalent to say that we take a perturbation of equation fl5^ which is analytic in 
a neighborhood of both r = and e = 0. System (l34l) provides an example where 
the upper bound of k = [(m — 1)/2J limit cycles bifurcating from the origin po 
of system (fTOj) . when only these perturbations are considered, is attained. In this 
way we have that Cycl(;fi"'^""^l,po) = [{m - 1)/2J. | 

Proof of Corollary I12L 

Let us consider system ffTOj) with a nilpotent focus at the origin and the correspond- 
ing Andreev number n. Let Vo{x, y) be an inverse integrating factor of this system 
which is analytic in a neighborhood of the origin. Therefore, we have an inverse 
integrating factor Vq*(x,?/) of the corresponding system f[T^ which is analytic in a 
neighborhood of the origin. We consider the Taylor development around r = of 
the following function: 

V*{rCse, r"Sn0) = </(^)r^^ + 0(r^+^), 

where vlj{6) ^ 0. 

Let us consider the transformation of system f|T3l) to an equation over a cylinder 
by generalized polar coordinates, see Lemma [311 and the corresponding inverse 
integrating factor V{r, 9) which is smooth and non-flat in r in a neighborhood of 
r = 0. We take the Taylor development around r = of the following functions: 

V{r,e) = v^{e)r^ + 0(r"^+i), Q{r,e) = 0„-i(^^)r"-i + 0(r"), 
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where Q{r,9) is the function defined in (^11) . We recall, see again the proof of 
Lemma EH that 9„_i(6') 7^ for any 9 e [0,T„). Indeed, by Lemma [201 we have 
that Vm{0) 7^ for any value of 6' G [0, T„). 

By statement (i) in Lemma [361 we deduce that vl,j{9) = Vm{0) 6„_i(6'). There- 
fore, we have that v*j^{6) 7^ for any value of 6* G [0,T„). On the other hand, 
since Vq{x, y) is analytic in a neighborhood of the origin, we deduce that v\j{6) is 
a (1, n)-quasihomogeneous trigonometric polynomial of weighted degree M. Using 
the symmetries of the generalized trigonometric functions described in statement 
(e) of Proposition [8], we see that 

4/(y) = (-I)^'-m(O). 

Therefore, we conclude that M needs to be an even number. 

We consider the value of m defined in Theorem [TOl and Remark [TTl and we have 
that M = 2n — 1 + m. Since M is even, we deduce that m is odd. Moreover, since 
the origin of system ffTOj) is assumed to be a focus, we have that m is an integer 
number with m > 1 and that m and n need to have the same parity, see statement 
(ii) of Theorem [TUl Thus, n is odd. | 

Proof of Corollary I19L 

As we have already stated, see Lemmas [25] and [31], in a neighborhood of these 

singular points, we can transform the system to an equation over a cylinder by 

means of (generalized) polar coordinates. Since the origin is a focus, we have that 

the periodic orbit r = is a limit cycle for the equation on the cylinder. Lemma 

[21] ensures the existence of a smooth, and non-fiat in r in a neighborhood of r = 0, 

inverse integrating factor V{r^ 9) for the equation over the cylinder. Indeed, we have 

that its vanishing multiplicity m at the origin is at least 1, that is, there exists a 

smooth function /(r, 9) defined on the cylinder such that V{r, 9) = r"* /(r, 9), with 

m > 1. The inverse integrating factor V{r,9) gives rise to an inverse integrating 

factor Vo{x,y) in cartesian coordinates. 

In the case of polar coordinates, we have that r = a/x^ + and 9 = arctan(y/x). 

We observe that the functions i^, and rf^ are bounded functions in a neigh- 

ax ' ay' ox ay ° 

borhood of the origin. By Remark [2] and if we consider system ([5]) with d > 1, we 
have that Vo{x,y) = r^^'^ f{r,9) with r and 9 expressed in cartesian coordinates 
and /(r, 9) a bounded function, with bounded derivatives, in a neighborhood of 
the origin. 

In the case of generalized polar coordinates, we have that r = x^" -f ny'^ 
and we observe that the functions jf- and r"~^|^ are bounded in a neighborhood 

ax chi CD 
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of the origin because 

If, in these expressions, we consider again the change of coordinates (fTS!) we have 
that 1^ = Cs^"~^6' and = Sn^. From the change of coordinates (ITS!) and 

using the definition of the generahzed trigonometric functions, see ([SD, it can be 
shown that 

d9 —ny d 

Thus, we have that the functions r|| and r'^^ are bounded in a neighborhood of 
the origin by an analogous argument as before. By Remark [TT| and if we consider 
system (fTOj) with n > 1, we have that Vo(x, ?/) = r'^+'^"--^ f(^r,6) with r and 6 
expressed in cartesian coordinates and /(r, 9) a bounded function, with bounded 
derivatives, in a neighborhood of the origin. 

We have, in both cases, that Vo{x,y) = r"' f{r,6), with r and 6 expressed in 
cartesian coordinates and /(r, 9) a bounded function, with bounded derivatives, in 
a neighborhood of the origin. Since m > 1, we have that the exponent a > 1 in the 
case of polar coordinates and a > n in the case of generalized polar coordinates. 
Thus, the limit of this product when (x, y) tends to the origin exists and it is 
equal to zero. Therefore, Vo{x, y) is continuous in a neighborhood of the origin and 
^o(0,0) = 0. 

By the chain rule, we have that 

dV^ ^ dV^dy^ ^ dVode_ ^ / a-ij. 0. ^ ^adl\ dr_ ^ d9_ 
dx dr dx dO dx \ ' dr j dx dO dx 

This expression can also be written as the product of a function tending to zero 
(because a > 1 or a > n, respectively, in each case) and a bounded function. Thus, 
the function ^ is continuous in a neighborhood of the origin and it is zero on 
the point (0,0). An analogous argument holds for We have that the function 
Vo(x, y) and its first derivatives are continuous and vanish at the origin. Therefore, 
Vo(x, y) is at least of class in a neighborhood of the origin. | 

We remark that in Example [TH we have given an inverse integrating factor 
Vo(x, y) for system (|T71) which might be only of class in a neighborhood of the 
origin, depending on the values of m and n. 
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